
figure;         
f = gcf;

colormap bone;
fontname(f,'Arial');

[X,Y] = meshgrid(alphavec,betavec);
contourf(X,Y,rmat,20,'LineStyle','none'); 
hold on
plot(alphax,betax,'p','MarkerFaceColor','white','MarkerEdgeColor','yellow','MarkerSize',13);
hold off

xlabel('Agglomeration in productivity (alpha)','fontsize',14);
ylabel('Agglomeration in amenities (beta)','fontsize',14);

ax = gca;        % define ax before using it
ax.FontSize = 14;

axis square

c = colorbar;
c.Label.String   = 'Relative Population of Location 1 to 2';
c.Label.FontSize = 14;
